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[i] We present an analysis of the seasonally varying ventilation of the stratosphere using 
one-way flux distributions. Robust transport diagnostics are computed using GEOSCCM 
subject to fixed present-day climate forcings. From the one-way flux, we determine the 
mass of the stratosphere that is in transit since entry through the tropical tropopause to its 
exit back into the troposphere, partitioned according to stratospheric residence time and 
exit location. The seasonalities of all diagnostics are quantified with respect to the month 
of year (a) when air enters the stratosphere, (b) when the mass of the stratosphere is 
partitioned, and (c) when air exits back into the troposphere. We find that the return flux, 
within 3 months since entry, depends strongly on when entry occurred: (34±10)% more 
of the air entering the stratosphere in July leaves poleward of 45°N compared to air that 
enters in January. The month of year when the air mass is partitioned is also found to be 
important: The stratosphere contains about six times more air of tropical origin during 
late summer and early fall that will leave poleward of 45° within 6 months since entering 
the stratosphere compared to during late winter to late spring. When the entire mass of the 
air that entered the stratosphere at the tropics regardless of its residence time is 
considered, we find that (51±1)% and (39±2)% will leave poleward of 10° in the 
Nothem Flemisphere (NF1) and Southern Flemisphere (SFT), respectively. 
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ventilation of the stratosphere: Robust diagnostics from one-way flux distributions, J. Geophys. Res. Atmos., 119, 293-306, 
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1. Introduction 

[ 2 ] The stratosphere-tropo sphere exchange (STE) of mass 
and tracers strongly influences the chemical and radiative 
properties of the upper troposphere and lower stratosphere 
[e.g., Morgenstern and Carver, 2001; Park et al., 2004; 
Jing et al., 2005; Hegglin et al., 2006; Pan et al., 2006]. 
Several studies have investigated the impact of STE on tro- 
pospheric species like ozone and nitrous oxide, with Liang 
et al. [2009] finding that as much as 78% of ozone in the 
Arctic upper troposphere is attributable directly to transport 
from the stratosphere. STE has also been shown to play a 
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key role in setting the springtime maximum of Northern 
Hemisphere (NH) ozone [Danielsen and Mohnen, 1977; 
Monks, 2000; Slohl et al., 2003] and in determining the 
seasonal cycles of tropospheric species including Strontium- 
90, Beryllium-7 [Fry et ah, 1960; Dibb et ah, 1992], 
chloroflourocarbons, and nitrous oxide [Nevison et al., 2004; 
Liang et ah, 2008]. 

[ 3 ] The tropical tropopause layer (TTL) is an important 
source of air for the stratosphere and provides the likely 
boundary condition for stratospheric water vapor [. Holton 
et ah, 1995]. In order to understand the key role that the 
TTL plays in setting the distributions of water vapor and 
other trace species, numerous studies have documented the 
“tape recorder”-like propagation of seasonally varying water 
vapor [e.g., Mote et ah, 1996; Schoeberl et ah, 2008b] and 
carbon monoxide [e.g., Rondel et ah, 2007; Schoeberl et 
ah, 2008a; Abalos et ah, 2012] from the tropical tropopause 
into the interior of the stratosphere. Few studies, however, 
have assessed the seasonality of stratospheric transport from 
the tropopause back to the tropopause. To the best of our 
knowledge, no studies have documented how seasonal trans- 
port variations at the tropical tropopause affect fundamental 
stratospheric transport time scales such as the time between 
entering and exiting the stratosphere. Transport through the 
TTL has a seasonally recurring component that is associated 
with the large-scale circulation and with convection; other 
variations occur more episodically. For example, volcanic 
injections of sulfur dioxide and aerosols into the stratosphere 
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occur infrequently but provide the largest natural perturba- 
tions of the stratosphere, with significant impact on global 
surface temperatures [McCormick et al, 1995]. In light of 
the key role that STE plays in impacting not only strato- 
spheric and tropospheric chemistry but also surface climate, 
it is important that its seasonality be quantified in a robust, 
reproducible manner. 

[ 4 ] Numerous studies have quantified STE in terms of the 
exchanges of ozone and water vapor because of their key 
roles as greenhouse gases and chemically active species in 
the upper troposphere and lower stratosphere [e.g., Olsen 
et al., 2003; Hsu et al., 2005; Tang and Prather, 2010; 
Tang et al, 2011]. The distribution of these species is 
the result of the complex interplay between transport and 
chemistry. Therefore, to isolate the role of transport, it is 
important to obtain a fundamental, quantitative understand- 
ing of STE using diagnostics that are tracer independent, 
that is, independent of the specific source/sink distribution 
of any particular chemical species. To this end, several stud- 
ies have focused on diagnosing cross-tropopause transport 
using one-way stratosphere to troposphere (S — > T) and tro- 
posphere to stratosphere (T -> S) air mass fluxes [e.g., St oh I 
et al, 2000; Wernli and Bourqui, 2002]. Hall and Holzer 
[2003] showed, however, that unless additional constraints 
are imposed on where air enters/exits the stratosphere and/or 
how long that air spends in the stratosphere (the strato- 
spheric residence time), these one-way fluxes are fundamen- 
tally ill defined because of an eddy-diffusive singularity at 
zero residence time. This singularity arises physically from 
the quasi-random nature of eddy diffusion at the shortest 
resolved length and time scales. Although it can be avoided 
by demanding that fluid elements stay for at least a thresh- 
old residence time on the other side of the tropopause before 
recrossing, the one-way cross-tropopause flux is in practice 
found to be a sensitive function of threshold time, rapidly 
increasing with decreasing threshold [Wernli and Bourqui, 
2002; Hall and Holzer, 2003]. In recognition of this, Orbe 
et al. [2012] introduced a robust numerical procedure for 
obtaining the one-way air mass flux partitioned with respect 
to stratospheric residence time, r, and with respect to the 
regions where air enters the stratosphere and returns to 
the troposphere. 

[5] This partitioned flux, or flux distribution, quantifies 
the ventilation of the stratosphere much along the lines of 
the work of Primeau and Holzer [2006] for the ventilation 
of the ocean. The one-way flux distribution thus quanti- 
fies the one-way flux without being rendered ill defined 
by the eddy-diffusive singularity. Note that the net air 
mass flux of Holton et al. [1995] and Appenzeller et al. 
[1996] is not affected by the singularity. However, because 
of eddy diffusion, large cross-tropopause tracer gradients 
can lead to significant tracer transport even if the net 
mass flux is zero, which is why the net mass flux is an 
incomplete diagnostic [e.g., Gettelman and Sobel, 2000; 
Nakamura, 2007]. 

[6] Here we extend the idealized study of Orbe et al. 
[2012] for statistically stationary flow to comprehensively 
modeled seasonally varying flow. Orbe et al. [2012] used 
integrations with a dry dynamical core run under perpetual 
Northern Hemisphere (NH) winter conditions to illustrate 
the latitude and residence time structure of the climatologi- 
cal zonal mean one-way flux across the thermal tropopause. 


Here our goal is to determine the one-way air mass fluxes 
and associated stratospheric residence times for a compre- 
hensively modeled realistic atmosphere, with a focus on 
seasonal variations. 

[ 7 ] We calculate one-way stratosphere-to-troposphere 
flux distributions in the Goddard Earth Observing System 
Chemistry Climate Model (GEOSCCM). The seasonal forc- 
ing, representation of moist thermodynamics, topography, 
and surface fluxes enable us to quantify STE for a realis- 
tic atmosphere in terms of the one-way air mass flux and 
diagnostics derived therefrom. Special attention is paid to 
quantifying the seasonality of the flux with respect to the 
month of year (a) when air enters the stratosphere, (b) when 
the mass of the stratosphere is partitioned, and (c) when air 
exits back to the troposphere. Our results provide a bench- 
mark for evaluating the tropopause-to-tropopause transport 
through the stratosphere in other models. 

2. Theory of One-Way Flux Distributions 

[ s | We quantify one-way STE across the thermal 
tropopause (see section 3.2). To track air parcels during 
their transit from where they enter the tropopause at region 
£2 ; , to the tropopause point of exit, r<^, we use the bound- 
ary propagator Green function Q (which has dimensions of 
inverse time) so that Q may be thought of as a tracer label 
that is applied to air on £2, at t, and removed at tq at time 
(tj, ti + d/,). Physically, Q(r, t|£2,, t i )dt l is the mass fraction of 
air at position r and time t that had last £2, contact during 
(tj, tj+dtj). The quantify Q also has the interpretation of being 
the age spectrum for air entering the stratosphere at the trop- 
ics at time t, [Hall and Plumb, 1994], and to the degree that 
the flow is statistically stationary, Q is also the distribution 
of transit times since last contact with the tropopause (i.e., 
transit time distribution, the distribution of air with times 
(t - ti, t-ti + d/,) since last contact with £2,). 

[ 9 ] The boundary propagator Q is computed as the pas- 
sive tracer response to a pulse in mixing ratio applied on 
£2, during (6, t, + d/,) that satisfies the source-free advection- 
diffusion equation 

|(pS) + V-J = 0, (1) 

where J is the advective-diffusive mass flux of Q and p is the 
density of air. For the generic case of advection with velocity 
u and Fickian diffusion with diffusivity tensor K, we have 
J = (pu - pKV)(/. The tropopause boundary conditions for 
Q [Holzer and Hall, 2000] are 

QOnJl^idi) = A(r Q ,£2 i)S(t-ti), (2) 

where 8(t - tj) is the usual Dirac delta function, and the sur- 
face mask A(rQ, £2 ; ) = 1 if lies on £2„ and A(i-q, £2,) = 0 
otherwise. Note that a boundary condition of zero applies to 
the entire tropopause for t > ti. This ensures that the Q label 
is removed on repeat contact with the tropopause. We will 
also refer to (/-labeled air as “£2, air”. 

[ 10 ] The one-way S — > T air-mass flux through ro that 
has resided in the stratosphere for a time in the interval 
(r, r + dr) is given by J(r^i,ti + r|£2„ t,)dr, where J is 
the component of J that is normal to the tropopause. Note, 
however, that we do not calculate these fluxes explicitly; 
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instead, we use a simple mask-budget approach [Orbe et al., 
2012] that is more robust and drastically simplifies our cal- 
culations (see experimental design, section 3). Thus, J is 
a flux-density distribution that partitions the flux of 52, air 
at t, and exiting at point rc 2 according to residence time r 
[Primeau and Hoher, 2006]. The boundary conditions (2) 
ensure that there is no return flux of 52, air from the tropo- 
sphere so that J is the one-way, or gross, S^T flux with 
dimensions of mass/area/time 2 . 

[n] While J is interesting in itself, Primeau and Hoher 
[2006] show that J can also be used to obtain the inte- 
grable distribution 'JZ defined so that 77dTd 2 ro is the mass 
of the stratosphere in transit at time t that entered at 52, 
and will exit at after a stratospheric residence time in 
the interval (r, r + dr). When 52, is the entire tropopause, 
In d 2 rQ / 0 °° d r 'JZ gives the entire mass of the strato- 
sphere at time t. The residence time distribution, IZ, (also 
“transport-mass distribution” [Holzer and Hall, 2008]) is 
given by 


1l(T,£li,rn;t)= f dti J(rnJi + (3) 

J t-T 

Physically, (3) states that the mass in transit at time t 
whose residence time will he in the interval (r, r + dr) con- 
sists of air that entered the stratosphere at times t-r through 
t and will leave the stratosphere a time r later at rQ . In the 
case of steady flow, (3) becomes IZ = r J , which states that 
if the steady flux of air whose residence time is in the inter- 
val (r + dr) acts for a time r, the mass of the stratosphere in 
that residence time interval will be flushed out [Primeau and 
Hoher, 2006]. Note that for annually repeating cyclostation- 
ary flow, 'JZ depends on t only through t = t mod 1 year, i.e., 
only on the time of year at which we partition the mass of 
the stratosphere. 

[ 12 ] Because IZ partitions the finite mass of the strato- 
sphere, it is integrable even when overlapping entry/exit 
patches introduce a singularity at r = 0. Therefore, IZ is use- 
fully summarized by its moments. The mass, /i, in transit 
at time t from 52, to regardless of residence time, r, is 
given by 


/ftA.r n;t) 


d 2 rn 



dr 7Z(t, Cli,rci; t), 


(4) 


where d 2 r<^ is the horizontal area of the grid box at r<^ and 
the mean 52, — >• residence time is given by 

T(S,-,rn;0= — — -d 2 r Q / dr r72fr,£2 !; r a ;t). (5) 

li(Qi,rnri) Jo 

Orbe et al. [2012] point out that f provides an upper 
bound on the time available for the mass of the strato- 
sphere in transit from 52, to rQ at time t to be subjected to 
stratospheric chemistry. 

[13] The flux distribution J is the fundamental transport 
diagnostic at the core of our analysis. Using a comprehen- 
sive general circulation model, we provide the first analysis 
of quantitatively realistic one-way STE in tenns of one-way 
flux distributions, the mass in transit from 52, — > ro, and the 
mean residence time of this mass. We focus, in particular, on 
the seasonal variations of these transport diagnostics. 


3. Experimental Design 

3.1. The Model 

[ 14 ] We use GEOSCCM Version 2, which couples the 
GEOS5-general circulation model [Rienecker et al., 2008] 
with a comprehensive stratospheric chemistry package 
[Douglass et al., 1996]. The model has a horizontal resolu- 
tion of 2° latitude by 2.5° longitude, with 72 vertical levels 
extending from the surface to 0.01 hPa. The model is forced 
with annually repeating 2010 greenhouse gases (GHG) and 
ozone depleting substances under the A1B and Al sce- 
narios, respectively [Intergovernmental Panel on Climate 
Change, 2001; World Meteorological Organization (WMO), 
2007]. Sea surface temperatures and sea ice concentrations 
are years 2000-2019 time averages taken from an integration 
of National Center for Atmospheric Research Community 
Climate System Model (CCSM) Model 3.0 subject to A1B 
GHG forcing. Solar forcings are held constant, and there is 
no representation of the quasi-biennial oscillation (QBO). 

[15] Spin-up to a statistically stationary state takes 10 
years after which we introduce our diagnostic tracers that 
are passively advected using a flux-form semi-Lagrangian 
scheme [Lin and Rood, 1996]. There is no explicit tracer dif- 
fusion term. Compared to other chemistry climate models, 
GEOSCCM performs well in terms of its mean stratospheric 
dynamical and thermal structure, trace gas distributions, 
and their decadal changes in the recent past [Eyring et 
al., 2010]. While the model underestimates tropical inter- 
annual variability due to the lack of a QBO, this is not 
expected to impact significantly on mean tropical transport 
[Li et al., 2012]. 

3.2. The Diagnostic Tracers 

[ 16 ] The thermal tropopause is defined using the standard 
World Meteorological Organization definition [WMO, 1957] 
and computed online for each time step. Because the main 
pathway from anywhere in the troposphere into the upper 
stratosphere is through the tropical tropopause, we con- 
sider only a single entry region, 52,, namely the tropopause 
from 10°S to 10°N, which is characterized by annual mean 
upwelling in the tropical lower stratosphere. The boundary 
conditions (2) are enforced by broadening the delta function 
to a square pulse of duration St, = 30 days and amplitude i. 
In practice, A(r<^, 52,), defined in (2), is extended for numer- 
ical convenience throughout the depth of the troposphere 
[Orbe et al., 2012]. In the case of multiple tropopauses, our 
procedure captures the lowermost tropopause. This suffices 
because our resolution is too coarse to resolve all tropopause 
folds. It is worth noting, however, that our diagnostics can 
be applied to any surface that separates two domains of the 
atmosphere, even if that surface is multiply connected with 
isolated bubbles of air from one domain embedded in the 
other domain. 

[ 17 ] Here we explore J for seasonally varying flow so 
that J depends explicitly on (a) when air enters the strato- 
sphere (i.e., entry time f,) and (b) when air exits the strato- 
sphere (i.e., exit time tj). Figure 1 shows a schematic of a 
fluid element traversing the stratosphere. This fluid element 
is in transit at time t at an interior stratospheric point r, has 
had a transit time (5 = t - 1, since last tropical tropopause 
contact, and leaves the stratosphere at tf after a residence 
time r = tf—tj. 
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Figure 1. Schematic illustration of tropopause-to-tropopause transport through the stratosphere. Air 
enters the stratosphere at time ti, is in transit in the interior of the stratosphere at time t and exits back 
into the troposphere at time tf. Air in the interior of the stratosphere at time t has a transit time £ = t - 1, 


since last contact with the tropopause; air leaving the 
r = tf- t t since entering at the tropopause. 

[is] In order to probe the seasonal cycle of the one-way 
flux with respect to entry time t„ we compute an ensem- 
ble consisting of three members. The first ensemble member 
consists of six pulse tracers, Q(r, f|£2,-, tf), where t, = January, 
March, ..., November in the first year of the model inte- 
gration following spin-up. The second and third members 
each consist of six such pulse tracers for the same months 


stratosphere has resided in the stratosphere a time 


in years 2 and 3 of the same integration. Following Li et al. 
[2012], who show that an ensemble of six pulses is sufficient 
to capture the seasonality of lower stratospheric age spectra 
in GEOSCCM, we assume that six pulses are also sufficient 
for capturing the seasonal variability of J . Our choice of 
the three ensemble members rests on the assumption that the 
largest variability is primarily seasonal and that interannual 


1 

m f ) 



d 2 r^i 

U f = [60°N, 90°N] 


[kgm 2 s 2 ] 



Figure 2. The ensemble mean one-way flux of air that entered the stratosphere through the tropical 
tropopause G, (a zonal strip between ±10°) at time t, and leaves poleward of 60°N, averaged over 
Qf = [60°N, 90°N]. (a) The seasonally varying one-way flux is decomposed into (b) a residence time- 
dependent annual mean component (averaged over entry time t, at fixed residence time r = tf- tf) and (c) 
deviations therefrom. 
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[ d 2 Ta f —G(T n ,ti+£\ni,ti) 

Jn t Jo 9 




Figure 3. (a) The ensemble-averaged boundary propaga- 
tor, Q, mass integrated over the latitudes poleward of 45° 
(NH solid, SH dashed) for air that entered the stratosphere 
in January (blue) and July (red), (b) The one-way flux dis- 
tribution, J , integrated over the tropopause poleward of 
45°N (NH solid, SH dashed). Shading in Figures 3a and 3b 
denotes one standard deviation from the ensemble mean. 

variability is small compared to seasonal variations. This is 
a reasonable assumption given the lack of a modeled QBO. 

[ 19 ] For each 30 day pulse, we compute the mass flux per 
unit residence time, J , as the mass of £2, air that crosses 
back into the troposphere per time step between succes- 
sive applications of the boundary condition (2) [Orbe et al, 
2012]. Equation (1) is integrated for 20 years (following the 
10 year spin-up to climatology before our diagnostic trac- 
ers are added), consistent with studies by Schoeberl et al. 
[2003] and Li et al. [2012], who show that about 20 years are 
needed to capture the long time scales of stratospheric trans- 
port. Integral (3) is evaluated after linearly interpolating J 
to monthly /, values for residence times greater than 1 month 
and to daily /, values for smaller residence times so that the 
small-r singular behavior is resolved and its contributions to 
[ i . and f accurately captured. The contributions to the inte- 
grals in (4) and (5) for r > 30 years account for roughly 1% 
of the total residence time integrated flux and are evaluated 
analytically by extrapolating the exponential tails of Q. 

3.3. Seasonality of J: Variations With Respect 
to Entry Time tj 

[ 20 ] In order to illustrate the nature of the seasonality of 
J, we show in Figure 2 the ensemble-averaged and zonally 
averaged one-way flux distribution, integrated over the NH 
high latitudes, in the (t,, tj) plane. The point of Figure 2 is 
to demonstrate that the annually averaged flux distribution 
is defined as an average over the annual cycle with respect 
to b at fixed residence time r = tf-tj. The flux denoted by 
J can then be decomposed into this mean (Figure 2b) and 


deviations therefrom (Figure 2c). Note that the deviations 
from the annual mean sum to zero along lines of constant r 
(the slanted contours of Figure 2b). 

[ 21 ] In Figure 2b, we introduce notation that designates 
annual means with respect to the time variable t x as A(t ) * = 
j f Q d t A(t), where t x is the time variable (either t, or tj) 
over which we are averaging at fixed r, with T= 1 year. 


4. Results 

[ 22 ] We present results derived from the ensemble- 
averaged boundary propagator. Recall that we refer to air 
entering the stratosphere at the tropics as “£2, air”. When 
emphasizing the dependence on entry time, we also refer to 
this air as “t, air” (e.g., “January air”). When we assign an 
uncertainty ±rr x to quantity X, cr x represents the standard 
deviation of X from its three -member ensemble mean. 

4.1. Evolution of Q 

[ 23 ] We begin by focusing on January and July air to 
explore the broad features of the seasonality of STE with 
respect to entry time tj. Any significant details that differ for 
other t, are noted below. 

[ 24 ] The evolution of these air masses is in accordance 
with our understanding of the stratospheric circulation: Air 
entering the stratosphere at £2, undergoes diabatic upwelling 
in the tropics and isentropic quasi-horizontal transport to 
high latitudes. We find that 2.9 x 10 18 kg more air enters the 
stratosphere in January than in July (33% more relative to 
the annual mean) in concert with strong seasonal upwelling 
at the tropical tropopause during January. 

[ 25 ] Figure 3 shows the evolution with transit time 
£ = t - 1, of the ensemble-averaged boundary propagator, 
Q( r, /|£2,, tj), mass integrated over the stratosphere poleward 
of 45°N and 45°S. These latitudes are chosen to highlight 
the transport of £2, air to high latitudes; equatorward of 
these latitudes for small transit times quasi-horizontal trans- 
port dominates in the subtropics and rapid eddy-diffusive 
transport dominates in the tropical entry region. 

[ 26 ] Figure 3a shows that (/-labeled tropical July (JUL) air 
peaks at a transit time of about 3 weeks in the NH, whereas 
January (JAN) air peaks at about 4 months. The fact that air 
that entered in July is more rapidly transported to NH high 
latitudes is consistent with the boreal summertime weak- 
ening of the subtropical mixing barrier that separates the 
tropics from the extratropics in the lower stratospheric over- 
world [Plumb, 1996; Neu and Plumb, 1996; Li et al, 2012]. 
Several studies have shown that this barrier is most perme- 
able during NH summer [Konopka et al., 2009; Bonisch et 
al, 2009], enabling air that is confined to the tropics to be 
spread poleward by eddies. This increased permeability has 
been associated with the extension of the Asian monsoon 
summer flow well into the lower stratosphere [Dunkerton, 
1995; Dethof et al, 1999; Haynes and Shuckburgh, 2000; 
Randel and Park, 2006]. 

[ 27 ] After peaking within 1 year of entering, £2, air decays 
exponentially with transit time. This decay quantifies the 
asymptotic return of stratospheric air to the troposphere 
where it is unlabeled [e.g., Prather, 1996; Ehhalt et al, 
2004; Pan et al, 2006; Holzer, 2009b]. Similar to the work 
in Li et al. [2012], we find that the decay e-folding time of 
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(a) 


[ d£ [ —G(Y n ,ti+£\ni,ti) (b) / dr + T\fl z ,ti) 

J J o 9 J n 


U = JAN 


U = JUL 


U = JAN 


U = JUL 



£ G [14 days, 3 months] 

[12 x 10 17 kg (11%)] [9.5 x 10 17 kg (13%)] 


r G [14 days, 3 months] 



£ G [3 months, 12 months] 


t G (3 months, 12 months] 

[1.8 x 10 10 kg s- 1 ] [9.0 x 10 9 kg s^ 1 ] 


Figure 4. (a) The ensemble-averaged and column-integrated boundary propagator Q, mass-averaged 
over air that has resided in the stratosphere (top) 2 weeks to 3 months and (bottom) 3 to 12 months, (b) 
The ensemble-averaged one-way S — > T air mass flux for air residing in the intervals shown. The globally 
integrated flux density integrated over residence time band r e [zi, r 2 ] is shown at the top of each map. 
Air that entered the stratosphere in January and in July are shown in Figures 4a and 4b, respectively. 
Numbers at the top of each map give the mass of air and the percentage of the incoming air mass that had 
last £2, contact in the transit time interval indicated. 


this mode is ~ 3 years. Figure 3 shows that within 3 years, 
Q has decayed to 25% of its initial value and variations with 
entry time t, have diminished substantially. Because trans- 
port into the SH is weaker than into the NH within the 
first 3 years of entering the stratosphere (Figure 3a), less air 
returns to the troposphere in the SH and interhemispheric 
differences rapidly even out during this time. 

[ 28 ] Figure 4 shows maps of the ensemble-averaged and 
column-integrated boundary propagator, integrated over var- 
ious transit time bands, to highlight geographic variations 
in the transport of 12, air (Figure 4a). Within the first 3 
months since entry into the stratosphere (Figures 4a, top) 
12, air is largely confined between the tropics and subtrop- 
ics, and meridional gradients are largest at the tropospheric 
jets, particularly in the NH where a greater mass of 12, air 
remains because cross-tropopause transport is strongly sup- 
pressed during NH winter [Chen, 1995]. Both JAN air and 
JUL air have prominent minima during NH summer and are 
centered on the Tibetan Plateau (Figure 4a, top right and bot- 
tom left). This feature coincides with the Asian monsoon and 
is consistent with observations showing localized extrema 
in upper troposphere/lower stratosphere tracer fields during 
NH summer, including anomalously high ozone [Randel and 
Park, 2006] and anomalously low water vapor [Rosenlof et 
al, 1997], methane [Park et al, 2004], and carbon monox- 


ide [Li et al, 2012]. The Asian monsoon is a deep circulation 
with closed anticyclones extending up to at least 70 hPa in 
July [Dunkerton, 1995]; at this level air remains confined 
for several weeks compared to only a fews days in the mid- 
dle troposphere and lower stratosphere [Randel and Park, 
2006]. This confinement suggests that £2, air trapped in the 
stratospheric monsoon flow has relatively ready access to 
the tropopause and hence remixes with the troposphere more 
quickly, resulting in the minima visible in Figure 4. 

4.2. The One-Way Flux Distribution, J 

[ 29 ] The fundamental quantity at the core of this 
study is the ensemble-averaged one-way flux distribution, 
J (rr 2 , //[ £2,, tj), which is simply the flux of Q across the ther- 
mal tropopause per unit residence time induced by boundary 
condition (2). In this section, we first discuss the seasonal- 
ity of the one-way flux with respect to entry time f„ and we 
then consider how the spatial distribution of J varies with 
exit time tj. 

[ 30 ] We focus first on the behavior of J for small res- 
idence times. The one-way flux of JAN and JUL air with 
residence times r — tf— t, of 3 months or less makes up 
(36±6)% and (55±3)% of the total flux poleward of ±45° 
(i.e., the return flux regardless of residence time) back into 
the troposphere. These small-r fluxes are larger in the NH 
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Figure 5. The ensemble-averaged and zonally averaged one-way S — *■ T mass flux , J tit,, that entered 
the stratosphere during a 30 day time interval tit, in (a) January and (b) July as a function of exit latitude 
and stratospheric residence time, r. Note the two color scales: (bottom) one for air with residence times 
less than 2 years and (top) the other for residence times greater than 2 years. Two linear color scales are 
used in order to capture the full dynamic range of J , including the large fluxes at small residence times. 
The residence time axis is stretched for r < 2 years so that the large short-r fluxes, and their rapid decay 
with r, are visible. The vertical lines mark the edges of f2,-. 


than in the SH, as shown in Figure 3b, where J has been 
separately integrated over each hemisphere poleward of 45°. 
Consistent with the fact that air is more rapidly transported 
to NH high latitudes in summer (see section 4.1), we find 
that the one-way flux across the NH tropopause for air with 
r < 3 months is (34 ±10)% larger for JUL air than for JAN 
air. The small-r fluxes back into the troposphere in the SH 
are larger for JAN air, as expected for SH summertime trans- 
port, but the amplitudes of J for both air masses are overall 
larger in the NH than in the SH so that, as in section 4.1, we 
mainly focus on the NH. 

[ 31 ] The maps of J in Figure 4b also reveal that the 
longitudinal structure of one-way STE depends strongly 
on how long the air mass has resided in the strato- 
sphere when r < 1 year. Air with r < 3 months is rapidly 
transported diffusively back across the tropical tropopause 
or quasi-horizontally across the tropopause of adjacent lat- 
itudes (Figure 4b, top). The more familiar storm track 
pattern, evident in the STE analyses of Seo and Bowman 
[2002], Sprenger and Wernli [2003], Tang and Prather 
[2010], and Tang et al. [2011], manifests itself for r >3 
months (Figure 4b, bottom). A few smaller-scale features 
are also worth pointing out: For JUL-air, enhanced one-way 
fluxes over Asia and Europe point to the key role that the 


Asian monsoon plays in ventilating the stratosphere during 
boreal summer. 

[ 32 ] Figure 5 shows the zonally averaged flux density 
J as a function of exit latitude and residence time. This 
figure shows that for a few weeks after first entering the 
stratosphere, J is overwhelmingly large at the entry region 
(±10°), as air readily diffuses back through the tropical 
tropopause, rendering the flux there difficult to estimate 
robustly unless diagnostics like J are used that partition 
the flux explicitly with respect to residence time. (Note that 
this diffusion is dominantly numerical here, because subgrid 
diffusion is not explicitly modeled.). A z 2 singularity at 
r = 0 is expected for Fickian diffusion in the continuum limit 
[Hall and Holzer, 2003; Primeau and Holzer, 2006; Orbe et 
al., 2012]. In the appendix, we show that the small-r power 
law is independent of entry time t,. With increasing r, these 
large fluxes decay rapidly, and within r ~ 200 days, J has 
decayed by 2 orders of magnitude. 

[ 33 ] For r > 2 years, away from the dominance of the 
nearly singular behavior, the pattern of J is strikingly dif- 
ferent at midlatitudes than in the tropics. The one-way flux 
returning across the midlatitude tropopause is made up of 
air that, at different stages during upwelling in the tropics, 
crosses the tropical barrier isentropically and spends a broad 
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Figure 6. The ensemble-averaged and zonally averaged one-way S — » T air mass flux for air entering 
the stratosphere at times /, and exiting back to the troposphere at times tj, averaged over ft, spanning (a) 
the SH high latitudes, (b) the SH midlatitudes, (c) the tropics, (d) the NH midlatitudes, and (e) the NH 
high latitudes. The vertical axis extends to 4 years in order to capture the strong seasonality of J with 
respect to t,. 


range of times recirculating in the extratropical stratosphere 
before descending back into the troposphere [Gettelman and 
Sobel, 2000; Sea and Bowman, 2002; Holzer, 2009a; Orbe 
et al., 2012]. Note that large short-r fluxes at midlatitudes 
are associated with rapid isentropic transport to midlati- 
tudes, not with the rapid eddy-diffusive near-singularity, 
because the midlatitudes do not overlap with the tropical 
entry region. At larger residence times r, extratropical fluxes 
decay exponentially to zero (not shown) as all of the 52, air 
eventually returns to the troposphere. The e-folding time for 
this decay is ~ 3 years. 

[ 34 ] To further quantify the seasonality of J , it is reveal- 
ing to examine the evolution of the one-way flux with 
respect to entry time t, and exit time tj. Figure 6 shows 
J(rn, f/|52,, t,) integrated over five key latitude bands. Large 
fluxes at midlatitudes are visible for r < 1 year in both 
hemispheres for late summer to early autumn entry times 
(Figures 6b and 6d). The one-way flux at high latitudes also 
demonstrates strong sensitivity to t, at small residence times 
because quasi-isentropic transport to polar regions is sup- 
pressed when a strong polar vortex is present. Specifically, 
polar fluxes are 1 order of magnitude smaller for winter 
entry times (Figures 6a and 6e) than for summer entry times 
(within each hemisphere). 

[ 35 ] While the flux of air with r < 2 years is seen in 
Figure 6 to depend strongly on the one-way flux of older 
air is only weakly dependent on f,. July maxima with respect 
to exit time tj are seasonally recurring in the tropics and may 


be explained in terms of the boreal summertime weaken- 
ing of the tropical-extratropical mixing barrier [Bonisch et 
al, 2009; Li et al, 2012]. (The seasonally recurring fluxes 
back into the tropics were noticeably absent in the idealized 
work of Orbe et al. [2012], because the flow was statis- 
tically stationary and representative only of perpetual NH 
winter.). In the subtropics, fluxes peak in late winter to early 
spring (Figures 6b and 6d) and at high latitudes during late 
summer (Figures 6a and 6e), independent of entry time. We 
explain this seasonality in exit time tf as follows: In late 
winter/early spring the presence of a strong polar vortex 
suppresses mixing to high latitudes, which leads to anoma- 
lously large fluxes out of the subtropics at the tropopause 
(i.e., maxima from late winter to early spring). At high lat- 
itudes during summer, when the polar vortex is weakest, 
tropical and subtropical air is free to move into high latitudes 
without being shunted to low and midlatitudes (i.e., max- 
ima from late summer to early autumn). Li et al. [2012] used 
similar reasoning to explain the seasonality of stratospheric 
age spectra. 

4.3. The Mass of the Stratosphere in Transit 
From ft, — *■ rfi and its Mean Residence Time 

[36] To further distill the rich information contained in the 
one-way flux distribution, we now focus on the partitioning 
of the mass of ft, air in transit during time of year 7 accord- 
ing to where the air eventually exits and its residence time 
in the stratosphere. From this transport-mass distribution we 
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Figure 7. The ensemble-averaged and zonally averaged mass of the stratosphere (per unit residence 
time) that is in transit from £2 , to rr> during month of year t and will leave the stratosphere after a residence 
time t, averaged over exit regions £2/ spanning (a) the SH high latitudes, (b) the SH midlatitudes, (c) the 
tropics, (d) the NH midlatitudes, and (e) the NH high latitudes. The vertical axis extends to 4 years in 
order to capture the strong seasonality of 1Z with respect to t. 


compute p, the mass in transit from £2, to tq regardless of 
residence time, and r, the mean residence time (the norm 
and first temporal moment). 

4.3.1. Transport-Mass Distribution, 1Z 

[ 37 ] Figure 7 shows the £2, — £2/ transport-mass dis- 
tribution 1Z (see equation (3)) as a function of the month 
of year t during which the £2, air mass is being partitioned 
and as a function of its eventual residence time, r. We use 
the same exit regions £2 f as for Figure 6, for comparison 
with J . Interestingly, we find that the £2, air mass des- 
tined to return through these £2/ regions depends strongly 
on t only for r <6 months. For t in late summer to early 
fall, air with r < 6 months (that is, air that entered the 
stratosphere anytime from late winter to late summer) is 
about 6 times more likely to be flushed back into the tro- 
posphere than for t in winter to early spring (i.e., air that 
first entered the stratosphere anytime from summer to late 
winter). Furthermore, this rapid ventilation is not hemispher- 
ically symmetric: Approximately 30% more mass leaves the 
NH than the SH within the first 6 months. Correspondingly, 
the amplitude of the seasonal cycle of 7 Z with respect to t is 
about two thirds smaller in the SH. For r > 6 months, how- 
ever, it does not matter during which month we partition the 
stratosphere: / Q d 2 rr 2 7?.(r, fi,, rc 2 ; t ) is largely independent 
of t for all exit regions £2/. 


4.3.2. Total £2, — ► rj 2 Mass in Transit, fi, and Mean 
Residence Time, r 

[38] To summarize the information contained in 1Z, we 
calculate /r(£2„ r q; ?), the mass of the stratosphere in transit 
from £2, to rr 2 during month of year t, regardless of resi- 
dence time defined by equation (4). Figure 8 (left) shows 
maps of the rf 2 dependence of //(£2,, rc 2 ) together with their 
zonal averages. It turns out that the 1 dependence of // is 
weak. We therefore focus on the annual average with respect 
to t because 7 Z itself depends strongly on t only for r < 
1 year. One key result evident from Figure 8 is that the 
largest mass of the stratosphere that enters in the tropics will 
idtimately leave primarily where (i) it entered (10)% and (ii) 
through isentropic pathways at midlatitudes, with (51±1)% 
and (39±2)% leaving between 10°N-90°N and 10°S-90°S, 
respectively. Except over £2„ where the large values of p 
are associated with the short-r diffusive singularity [Orbe 
et al., 2012], the largest masses return to the troposphere 
downstream of the storm tracks. 

[ 39 ] It is useful to contrast p with previous diagnos- 
tics of STE. While the S — > T mass fluxes diagnosed 
from Lagrangian trajectories by Seo and Bowman [2002] 
and Sprenger and Wernli [2003] share a broadly similar 
geographic pattern with p, their mass fluxes and // are fun- 
damentally different diagnostics: The fluxes calculated in 
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Figure 8. (left) Map of the exit location dependence of the ensemble-averaged mass, /x, of the strato- 
sphere in transit from Q, to exit location ro, independent of when that air is in transit in the stratosphere. 
We have plotted the mass leaving the stratosphere at tq per unit area (m 2 ) so that, when integrating over 
the Earth’s surface, one recovers the mass leaving the stratosphere that came in at £2,-. (right) As in the 
left, except now shown for the zonal mean and plotted in terms of the mass leaving per degree latitude. 


those studies are computed for trajectories that start in the 
interior of the stratosphere, while // captures the mass of the 
stratosphere in transit between successive tropopause cross- 
ings. It is also worth emphasizing that /x does not suffer from 
any singularities, and unlike the one-way flux regardless of 
residence time is well defined for all rr 2 . Mathematically, 
lZ{r, £2„ rc 2 ; t) is integrable with respect to r for all ro . 

[ 40 ] We now consider the mean residence time, f, of 
£2, — > ro air, defined by equation (5). The mean residence 
time turns out to have very little variation with longitude 
because the norm ji is divided out so that the mass whose 
residence time is computed does not come into play. We 
therefore consider only zonal averages of f, which are 


shown in Figure 9. As expected from the singular behav- 
ior of J and 1Z for small r and ro € £2„ Figure 9a 
shows that f <100 days over the tropical entry region, 
which reflects the overwhelming importance of the eddy- 
diffusive return flux at short residence times. (This behavior 
was also captured in the idealized experiments of Orhe et al. 
[2012]). Mean residence times increase with latitude to ~ 5 
years poleward of ± 30° and are about 1 year shorter in 
the NH than in the SH, reflecting the more vigorous flow 
and residual circulation of the NH. While it is tempting to 
attribute these shorter mean residence times to differences 
between the strength of the residual mean circulation of each 
hemisphere, we emphasize that this can only be a partial 



Figure 9. (a) The ensemble-averaged and zonally averaged mean residence time, f , of Q, air, annually 

averaged over the seasonal cycle with respect to t, the time of year when the air is being partitioned, 
(b) Tropopause variability quantified in terms of the standard deviation of tropopause pressure from the 
annual mean, (c) Ensemble -mean deviations of the seasonally varying mean residence time, [r], from the 
annual mean in Figure 9a. Contours are spaced 25 days apart. 
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explanation. The residual stream function is not sufficient 
to characterize stratospheric transport, a point that is now 
widely appreciated in the stratospheric dynamics community 
[Plumb, 2007], 

[ 41 ] Finally, Figures 9b and 9c show that the residence 
time of 52, air depends modestly on t. This seasonality with 
respect to when the air mass is partitioned is strongest for 
air leaving the stratosphere in the subtropics: The 52, air 
in transit in summer that is destined to leave in the sub- 
tropics will have resided in the stratosphere ~ 200 days 
longer than such air in transit during the other seasons. One 
can understand this result by considering integrating TZ with 
respect to r. The distribution TZ, shown in Figure 7, depends 
most strongly on t for r <1.5 years and for exit latitudes 
spanning ±(10°, 40°) (Figures 7b and 7d). Because mean 
residence time, r, preferentially weights the contributions 
made by older air masses, f is only weakly seasonal outside 
of the subtropics. Furthermore, the residence time distribu- 
tion for air that leaves the stratosphere at the subtropics is 
skewed to larger residence times for t in early spring to 
early summer compared to other times of the year. For that 
reason, the 52, air in transit in summer that is destined to 
leave in the subtropics will have resided in the stratosphere 
~ 200 days longer than air in transit during other seasons 
(Figures 9c). 

[ 42 ] It is important to distinguish mean residence time, 

r, from stratospheric mean age evaluated at the tropopause, 
r n , which is defined as the first moment of the distribution 
of transit times from 52, to repeat contact with the tropopause 
[Holzer et al, 2012]. While both time scales are average 
transit times between entry and exit at the tropopause, Holzer 
et al. [2012] show that T^ is an average over the population 
of fluid elements that are exiting the tropopause, whereas r 
is an average over all £2 — > fluid elements in the entire 

stratosphere. Because physically distinct populations of fluid 
elements are averaged for r q and r, the aggregated effects 
of eddy-diffusive mixing manifest themselves differently. 

5. Conclusions 

[ 43 ] We have quantified the seasonal ventilation of the 
stratosphere in terms of flux distributions that diagnose 
the transport of stratospheric air from entry at the tropi- 
cal tropopause to exit back into the troposphere. This has 
been accomplished with GEOSCCM, a comprehensive cir- 
culation model that provides a realistic representation of the 
dynamics and transport of the atmosphere | Byring et al., 
2010; Strahan et al., 2011]. The key quantity at the core of 
our analysis is the one-way flux distribution, J , that par- 
titions the cross-tropopause mass flux of air that entered 
the stratosphere between ± 10° according to where this 
air exits the stratosphere, tq, and according to its strato- 
spheric residence time, r. The one-way flux distribution has 
been shown to capture transport robustly without being ren- 
dered ill defined by the short- r eddy diffusive singularity 
[Hall and Holzer, 2003; Primeau and Holzer, 2006; Orbe et 
al, 2012], 

[ 44 ] Our previous work on one-way flux distributions 
exploited simplifications afforded by the statistical station- 
arity of the modeled flow [Orbe et al, 2012]. However, 
for realistic atmospheric flows, seasonal variations of the 
transport add considerable complexity. This seasonality is 


quantified here with respect to the month of year (a) when 
air enters the stratosphere, (b) when the mass of the strato- 
sphere is partitioned, and (c) when air exits back into the 
troposphere. Because of this added complexity, we consider 
only stratospheric air that entered in the tropics (£2, air), not 
the entire mass of the stratosphere as in Orbe et al. [2012]. 
From the one-way flux distribution, we compute the mass 
of the stratosphere undergoing 52, — > rc 2 transport and its 
stratospheric mean residence time. Our main findings are 
as follows: 

[ 45 ] 1. The evolution of £2, air during its first year in the 
stratosphere depends strongly on the time of year when entry 
occurs. For example, nearly twice as much JUL 52, air is 
transported to NH high latitudes during its first 3 months 
since entry compared to JAN £2, air. Correspondingly, north 
of 45°N the flux of 52, air back into the troposphere after a 
residence time r < 3 months is (34 ±10)% larger for JUL 
52, air compared to JAN 52, air. This is due to the fact that in 
the NH JAN 52, air tends to remain confined to low latitudes 
until the subtropical mixing barrier in the lower stratosphere 
weakens during boreal summer. The longitudinal structure 
of the flux for r < 1 year points to the key role that the Asian 
monsoon plays in rapidly ventilating the stratosphere during 
boreal summer. 

[ 46 ] 2. The flow rate of air from 52,, through the strato- 
sphere, and back into the SH troposphere is weaker than the 
flow rate back into the NH troposphere. However, because 
52, air within the SH is also slower to return to the tropo- 
sphere, the stratospheric mass of 52, air is distributed about 
evenly between the hemispheres 3 years after entry. After 3 
years, 75% of the 52, air has returned to the troposphere and 
its subsequent evolution is dominated by exponential decay 
with an e-folding time of ~ 3 years as it continuously leaves 
the stratosphere. The corresponding decay of the one-way 
return flux is modulated by a seasonally recurring compo- 
nent that peaks during summer in the tropics and during 
spring in the subtropics, moving to high latitudes by late 
summer. This seasonal phasing between the subtropical and 
high-latitude peak fluxes is consistent with the seasonal evo- 
lution of the polar vortex, which is strongest in midwinter 
when transport to high latitudes is suppressed. 

[ 47 ] 3. The 52, air mass regardless of residence time 
returns to the troposphere primarily at the midlatitude storm 
tracks (and also in the tropics due to near singular short- 
r diffusive fluxes). Averaged over 1, the time of year when 
the air mass in transit is partitioned, we find that (51±1)% 
of the 52, air mass is destined to return poleward of 10°N, 
and (39±2)% is destined to return poleward of 10°S, with 
the remainder returning in the tropics. The 52, air mass leav- 
ing in the tropics does so via the near-singular diffusive 
fluxes, which are purely numerical in our model and hence 
an aspect of the transport that we cannot reliably quantify. 
However, the qualitative aspects of these singular fluxes do 
match theoretical expectations (see Figure la), and they are 
consistent with modeling results in other contexts [Primeau 
and Holzer, 2006; Holzer, 2009a, 2009b; Orbe et al, 2012]. 

[ 48 ] 4. In contrast to the 52, air mass regardless of resi- 
dence time, the portion of this mass with residence times 
less than 6 months depends strongly on the time of year, t, 
during which the air mass is partitioned. Specifically, dur- 
ing September 6 times more 52, air with r = 6 months will 
leave the NH stratosphere compared to during March. For 52, 


303 


ORBE ET AL.: SEASONAL VENTILATION OF THE STRATOSPHERE 



Figure Al. Power law scaling of the area-integrated return 
flux J back through the tropical entry region between 10°S 
and 10°N for small stratospheric residence time, r. Six 
lines are shown for air entering the stratosphere in January, 
March, ..., November. The curves have been normalized by 
Jo, the value of J for r = 1 day. 

air with residence times longer than 6 months, however, this 
seasonality with respect to t is nearly negligible: During any 
month of year the stratosphere contains similar £2, air mass 
fractions with r > 6 months that are destined for a given 
exit region. Consequently, the £2, — »■ rQ air mass regard- 
less of residence time, //, has only weak seasonality with 
respect to t. 

[ 49 ] 5. The annually and zonally averaged mean residence 

time, f, of air undergoing £2, — > transport is roughly 

zero in the tropics and increases sharply with latitude to a 
nearly constant value of 5.7 ±0.2 years poleward of 45°S and 
5.1 ±0.1 years poleward of 45°N. Because f is independent 
of the mass for which it represents an average, variations of 
f with longitude are negligible. The residence time of £2, 
air has modest seasonality with respect to the month of year 
during which the air mass is partitioned. This seasonality is 
strongest for air leaving the stratosphere in the subtropics 
with an amplitude of ~ 20% of the annual mean. 

[ 50 ] The main point of this paper has been to apply a new 
tracer independent diagnostic of stratospheric transport to 
the comprehensively modeled flow of GEOSCCM. The flux 
distributions represent the one-way flux of air, and not of any 
particular trace species. However, to model specific trace 
gases, the flux distributions and their underlying boundary 
propagators can be convolved with tracer boundary condi- 
tions at the tropopause, weighted with decay functions of 
the time spent in the stratosphere. In this way our diagnos- 
tics can be used to disentangle the role of transport from that 
of chemistry. 

[ 51 ] Finally, recent studies have shown that future 
changes in STE will have important implications for the 
distribution of ozone [. Hegglin and Shepherd, 2009; Zeng 
et al., 2010], for the oxidizing capacity of the troposphere 
[e.g., Kentarchos and Roelofs, 2003], and for tropospheric 
air quality [e.g., Stohl et al., 2000; Cooper et al, 2005], It 
is therefore important to quantify how STE will change with 
the climate, and this will be explored in future work. 


Appendix A 

[ 52 ] Here we show the short-r dependence of 
r | £2 , , t : ) on r for £2/ = £2, for six different entry 

times t t into the stratosphere. Plotted is the area average of 
r|£2 j,tj) over £2/ = £2, = [10°S, 10°N] as a function 
of r for 1 day < r < 40 days. For clarity, all curves have 
been normalized by Jo, the flux at r = 1 day. The dashed 
line indicates a r 5 power law to guide the eye and shows 
that, remarkably, J diverges approximately like r 3/2 as 
r — > 0, as expected for Fickian diffusion [Hall and Holzer, 
2003; Orbe et al., 2012]. This implies that the combina- 
tion of resolved eddy diffusion and numerical diffusion in 
GEOSCCM results in simple downgradient diffusion. This 
“Fickian divergence” is seen to be insensitive to the time of 
year (7,) during which the air entered through £2,. 

[53] Acknowledgment. This work was supported by an NSF grant 
ATM-08547 1 1 (M.H. and L.P.). 
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